home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SPODI.z / SPODI
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSPPPPOOOODDDDIIII((((3333FFFF))))                                                            SSSSPPPPOOOODDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SPODI   - SPODI computes the determinant and inverse of a certain real
  10.      symmetric positive definite matrix (see below) using the factors computed
  11.      by SPOCO, SPOFA or SQRDC.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE SPODI(A,LDA,N,DET,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAA REAL(LDA, N)
  21.         the output  A  from SPOCO or SPOFA
  22.         or the output  X  from SQRDC.
  23.  
  24.      LLLLDDDDAAAA INTEGER
  25.         the leading dimension of the array  A .
  26.  
  27.      NNNN INTEGER
  28.         the order of the matrix  A .
  29.  
  30.      JJJJOOOOBBBB INTEGER
  31.         = 11   both determinant and inverse.
  32.         = 01   inverse only.
  33.         = 10   determinant only.  On Return
  34.  
  35.      AAAA If SPOCO or SPOFA was used to factor  A , then
  36.         SPODI produces the upper half of INVERSE(A) .
  37.         If SQRDC was used to decompose  X , then
  38.         SPODI produces the upper half of INVERSE(TRANS(X)*X),
  39.         where TRANS(X) is the transpose.
  40.         Elements of  A  below the diagonal are unchanged.
  41.         If the units digit of JOB is zero,  A  is unchanged.
  42.  
  43.      DDDDEEEETTTT REAL(2)
  44.         determinant of  A  or of  TRANS(X)*X  if requested.
  45.         Otherwise not referenced.
  46.         Determinant = DET(1) * 10.0**DET(2)
  47.         with  1.0 .LE. DET(1) .LT. 10.0
  48.         or  DET(1) .EQ. 0.0 .  Error Condition
  49.  
  50.      AAAA division by zero will occur if the input factor contains a zero on the
  51.      diagonal and the inverse is requested.  It will not occur if the
  52.      subroutines are called correctly and if SPOCO or SPOFA has set INFO .EQ.
  53.      0 .  LINPACK.  This version dated 08/14/78 .  Cleve Moler, University of
  54.      New Mexico, Argonne National Lab.  Subroutines and Functions BLAS
  55.      SAXPY,SSCAL Fortran MOD
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.